1 


NASA  Contractor  Report  177958 


CHEHOVISCOSITY  MODELING  FOR  THERMOSETTING  RESINS  -  11 


T.  H.  Hou 

PRC  Kentron,  Inc. 

Aerospace  Technologies  Division 
Hampton,  Virginia  23666 


(KASA-CR-177958)  CHPIMOVISCOSITY  MODELING  N85-352«3 

FOr  THFFMOS*!TTING  RPSISS,  2  Final  Report 
(PF.C  Kentron,  Inc.)  36  p  HC  NO  3/M  F  A01 

CSCL  11L  Onclas 

«3/2U  22295 


Contract  NASl-18000 


August  1985 


njASA 

Nat'ooal  Aeronautics  and 
Space  AdmiPrst  fatten 

Langtey  Research  Center 

Hampton,  Virginia  23665 


h:  'is 


;/ 


QUALITY  INSPECTED 


8 


19951128  100 


DEPa\RTMENT  of  defense 

PLASTlCb  TLCMt'>;iCAL  EVALUATION  CENTER 
AR D  cb"v; , '  "doTl K r  n7 ’  07501 


/ 

/ 


THIS  DOCUMENT  IS  BEST 
QUALITY  AVAILABLE.  THE  COPY 
FURNISHED  TO  DTIC  CONTAINED 
A  SIGNIFICANT  NUMBER  OF  _ 
PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


I 


i twi VAW V  ,• 


--  9  -  OESCRIPTIf  NOTE:  FINAL  REPT., 

--10  -  PERSONAL ilTHOfiS:  SAXTON  ,J.  C.  IKRAHER  ,N.  P.  :R0BERTS0N,D.  L.  ; 

FORTUNE,,!.  A.  :LEGGETT,N.  E.  ; 

"11  -  REPORT  DATE:  AUG  3C.  1974 

”12  *  pagination:  160P 
”14  -  REPORT  NUMBER: 


”  1  OF  1 

.  «4DTIC  does  NOT  HAVE  THIS  ITEHm 

i  ”  1  -  AD  NUMBER:  D439930 

;  ”  5  -  CORPORATE  AUTHOR:  NATIONAL  AERONAUTICS  AND  SPACE  ADMINISTRATION 
HAMPTON  VA  LANGLEY  RESEARCH  CENTER 

”  6  -  UNCLASSIFIED  TITLE:  CHEMOVISCOSITY  MODELING  FOR  THERMOSETTING 

RESINS-2.  , 

.  -  9-  DESCRIPTIVE  NOTE:  FINAL  REPT., 

”10  -  PERSONAL  AUTHORS:  HOU.T.  H.  ; 

.  ”11  -  REPORT  DATE:  AUG  ,  1985 

”12  -  PAGINATION:  36P 

•  ”15  -  CONTRACT  NUMBER:  NASl-18000 

”18  -  MONITOR  ACRONYM:  NASA 

-  ”19  MONITOR  SERIES:  CR-177958 

”20  -  REPORT  CLASSIFICATION:  UNCLASSIFIED 

”22  -  LIMITATIONS  (ALPHA):  APPROVED  FOR  PUBLIC  RELEASE:  DISTRIBUTION 
'  ”  UNLIMITED.  /WAILADILITY I  N.^HONA^-ri^HNIC.AL  INFOPHATIOli  C'gW‘'IO[, 

I  ”  <PRINtFlEkPi  VAi  OClul.  Ni:  iOiiC. 
f  ”33  -  LIMITATION  CODES:  1 

L  ■  , 

-END  .„.y  FOR  NEXT  ACCESSION  END 


Alt-Z  FOR  HELP3  ANSI  3  HDX  3  3  LOG  CLOSED  3  PRINT  OFF  3  PARITY 


"ir  «.»■' 


Chernov iscosity  Modeling  for  Theraosetting  Resins  -  II 

I  ' 


TABLE  OF  CONTENTS 

I  .  ' 

-  ‘  ABSTRACT  . .  . . 

FOREWORD  ■)  . . . .  .  .  .  . 

1*  -INTRODUCTION  ,  .  .  .  .  .  .  .  .  ......  .  .  , 

.  .  .  .  ....  ..... 

II.  THEORY  .  .  . 

III.  EXPERIMENTAL  DATA  . . . 

•••••••  *••••  •  •••4 

Degree  of  Cure  a(t,T)  ...  . 

Chernov iscosity  n(t,T)  . .  . . 

IV.  MODEL  PREDICTIONS  AND  DISCUSSION  . 

Isothermal  Case  ......... 

Dynamic  Heating  Case  . 

V.  CONCLUSIONS  ........... 

REFERENCES  .  . 

APPENDIX  -  GENERALIZED  NEWTON'S  METHOD  .  .  .  .  . 

LIST  OF  FIGURES  ......  . 

LIST  OF  COMPOSITES  PROCESSING  LABORATORY  RESEARCH  REPORTS  .  .  . 


PAGE 

ill 

V 

T 

3 

4 
7 
7 
9 

11 

11 

13 

18 

19 

20 
21 
22 


I  . 


ABSTRACT 

A  new  analytical  model  for  simulating  chemoviscosity  of  thermosetting 
resin  has  been  formulated.  The  model  is  developed  by  modifying  the 
well-established  Williams-Landel -Ferry  (WLF)  theory  in  polymer  rheology  for 
thermoplastic  materials.  By  assuming  a  linear  relationship between  the  glass 
transition  temperature  Tn(t)  and  the  degree  of  cure  a(t)  of  the  resin  system 
under  cure,  the  WLF  theory  can  be  modified  to  account  for  the  factor  of 
reaction  time.  Temperature-dependent  functions  of  the  modified  WLF  theory 
constants  Ci(T)  and  C2(T)  were  determined  from  the  isothermal  cure  data  of 
Lee,  Loos,  and  Springer  for  the  Hercules  3501-6  resin  system.  Theoretical 
predictions  of  the  model  for  the  resin  under  dynamic  heating  cure  cycles  were 
shown  to  compare  favorably  with  the  experimental  data  reported  by  Carpenter. 
This  work  represents  progress  toward  establishing  a  chemoviscosity  model  which 
is  capable  of  not  only  describing  viscosity  profiles  accurately  under  various 
cure  cycles,  but  also  correlating  viscosity  data  to  the  changes  of  physical 
p.-operties  associated  with  the  structural  transformations  of  the  thermosetting 
resin  systems  during  cure. 
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CHEMOVISCOSITY  MODELING  FOR  THERMOSEniNG  RESINS  -  II 
1.  INTRODUCTION 

There  are  two  different  basic  approaches  in  chemoviscosity  modeling  of 
ttiennosetttog  resin?  One  approecK,  typically  represented  by  the  « 

(toller  [1],  Is  to  empirically  formolate  a  rate  equation  uhich  relates  the 
Change  of  chemoviscosity  »ith  reaction  time.  For  a  resin  system  cured  under  a 
dynamic  heating  cure  cycle  (as  commonly  encountered  in  autoclave  processing 
for  composite  materials),  such  an  approach  usually  yields  a  model  which  is 
inadequate  in  describing  accurately  the  nonlinearity  of  chemoviscosity  as  a 
function  of  reaction  time.  As  Tajima  and  Crozier  [2]  pointed  out, 
modeling  approach  renders  itself  to  the  limitation  of  batch-specific.  The 
model  parameters  cannot  be  readily  related  to  the  chemical  and  rheological 
properties  of  the  reacting  systen  as  well. 

The  second  approach  is  based  upon  a  modification  to  the  well-established 
viscosity-temperature  relationship  existing  in  polymer  rheology  for 
thermoplastic  materials.  The  parameters  in  such  an  equation  can  be  expressed 
in  terms  of  polymerization  kinetics,  and  the  chemoviscosity  profiles  as  a 
function  of  reaction  time  can  then  be  modeled  for  thermosetting  resin  systems. 

The  applicability  of  the  modified  Williams-Landel -Ferry  (WLF)  theory  [3] 
in  chemoviscosity  modeling  for  thermosetting  resin  has  been  studied  by  Tajima 
and  Crozier  [2.4],  Apicella  et  al.  [5],  and  Hcj  [6]  among  others.  It  has  been 
extensively  documented  in  the  literature  that  temperature-dependent  viscosity 
of  thermoplastic  material  can  be  accurately  described  by  the  WLF  equations 


For  a 


within  100°K  above  its  glass  trat.‘ition  temperature  (Tq). 

,  \ 

thermosetting  resin  system  under  cure,  it  is  reasonable  to  assume  that  Tg(t) 
of  the  material  is  always  within  100“K  lower  than  the  cure  temperature.  The 
WLF  equation  should  then  be  applicable  to  every  instant  state  of  the  resin 
system  under  curing.  The  central  theme  in  modifying  the  WLF  theory  to 
describe  the  chemoviscosity  of  thermosetting  r^sin  is  to  take  into  account  the 
reaction  time  factor. 


In  an  earlier  paper  [6],  two  assumptions  were  proposed,  namely,  that  the 
rate  constant  of  reaction,  kj,  at  any  temperature  T,  is  (i)  diffusion 
controlled  and  is,  therefore,  inversely  proportional  to  the  viscosity 
nj(t)  of  the  reactive  medium,  and  (ii)  directly  proportional  to  the  rate 
of  change  of  glass  transition  temperature  Tg(t).  The  modified  WLF  equation 
becomes  a  first  order  ordinary  nonlinear  differential  equation.  Numerical 
solutions  have  also  been  shown  to  compare  favorably  with  the  experimental 
results  for  several  thermosetting  systems  under  isothermal  and  dynamic  heating 
cure  conditions.  It  has  been  concluded  that  the  flexibility  demonstrated  by 
the  modified  WLF  theory  [6]  can  be  conveniently  exploited  to  establish  an 
analytical  model  with  high  degree  of  accuracy  for  the  chemoviscosity  of  any 
thermosetting  resin  system  under  various  cure  cycles.  The  physical 
significances  of  the  material  parameters  selected  for  the  model  were,  however, 
difficult  to  extract  for  the  particular  resin  system  under  investigation. 


As  a  next  step  toward  the  improvement  of  the  chemoviscosity  modeling  for 
a  thermosetting  resin  system,  a  new  modification  of  the  WLF  theory  is  proposed 
in  this  report.  Comparisons  between  the  theoretical  predictions  and 
experimental  data  for  chemoviscosity  of  a  resin  system  (Hercules  3501-6)  under 
both  isothermal  and  dynamic  heating  cure  conditions  will  be  discussed. 


List  of  Symbols 


Aj,A2»A3 

a,b,c,d 

B 

C,(T).  C2{T) 


DSC 

AE  1  •  /£  2  *  3 

Kj,K2,K3 

R 

T 

Tc 

Tg{t) 

Tg{») 

^2 


t 

tf 

At 

WLF 

nj(  t ) 
a(t) 


Material  parameters  (min”^)  specified  in  Table  1 
Constants  used  in  Eq.  (5) 

Material  parameter  specified  in  Table  1 

Temperature-dependent  WLF  constants  defined  by  Eqs.  (7)  and 
'1(8)  " 

Differential  Scanning  Calorimeter 
Activation  energies  (J/mol)  specified  in  Table  1 
Material  constants  (min“M  defined  in  Eqs.  (3)  and  (4) 
Universal  gas  constant  (=1.987  Cal/gm-mole/®K) 

Temperature  (“K) 

Cure  temperature  ("K),  T^  =  Tg(co) 

Glass  transition  temperature  of  the  resin  system 
under  cure  ("K) 

Glass  transition  temperature  of  the  resins  system  at 
completion  of  cure  (“K) 

Glass  transition  temperature  of  the  uncured  resin  g 
system  (“K) 

Cure  time  (minutes) 

Time  (minutes)  defined  in  Eq.  (6) 

Time  required  for  completion  of  cure  (minutes) 

Time  interval  (minutes)  used  in  solving  Eq.  (9) 

Williams-Landel -Ferry 

Chemoviscosity  at  temperature  T  (poises) 

Degree  of  cure  of  the  thermosetting  resin  during  cure 
as  determined  by  DSC  thermal  analysis 


II.  THEORY 


The  William-Landel -Ferry  (WLF)  theory  [3]  is  given  by 


•Vr  ^  -Ci(T.T  ) 

-v;  9 


(1) 


for  Tg  <  T  <  Tg  +  100°K. 


and  Cj  are  two  material  constants,  and  ny  represents  viscosity  at 
temperature  T  of  the  given  thermoplastic  material  which  possesses  a  glass 
transition  temperature  Tg,  The  normal  use  of  the  WLF  equation  for 
thermoplastic  materials  requires  that  Tg  be  constant  while  the  temperature  T 
is  varied  for  the  specific  polymer  under  study.  However,  during  cure  of 
thermosetting  resins,  the  monomers  are  initially  polymerized  and  crosslinks 
are  formed  later.  This  is  a  system  where  Tg(t)  is  changing  and  the  curing 
temperature  T  is  held  constant  (in  an  isothermal  cure  case,  for  example).  The 
glass  transition  temperature  Tg  rises  continuously  and  may  eventually 
approach  the  curing  temperature.  Over  the  entire  curing  cycle,  the  material 
structure  actually  undergoes  continuous  phase  transformations  from  a  low 
molecular  weight  liquid  to  a  high  molecular  weight  polymeric  melt,  and 
eventually  transforms  to  a  crosslinked  network.  It  is  reasonable  to  assume 
that  Tg{t)  of  the  material  is  always  lower  than  the  cure  temperature  T,  and 
that  (T-Tg(t))  is  always  within  lOO^K.  The  WLF  theory  (Eq.  {!))  should  then 
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be  applicable  to  all  of  the  different  polymeric  structure  phases  during  cure. 
Before  Eq.  (1)  canibe  applied  to  describe  the  chamoviscosity  of  the  curing 
resin,  however,  modifications  need  to  be  made  in  order  that  reaction  time,  t, 
of  the  thermosettl/ng  resin  system  1s  properly  taken  Into  account.  This  can  be 

accomplished  as  follows: 

■I  '  ’  ■ 

"  ■  ^  that  Tg(t)  is  a  linear 'function  of  the  degree  of 

cure  ?(t)  of  the  thermosetting  resin  can  be  made,  which  o’ves: 

T,(t)  •  Tj  .  (T^-TJ)  .(t)  (2> 

where  T®  =  Tg{0)  is  the  glass  transition  temperature  of  the  unreacted  resin 
system,  Is  the  cure  temperature  In  “K,  and  t  Is  the  reaction  time  in 
minutes.  Equation  (2)  should  satisfy  the  boundary  conditions,  namely, 

T  =T®,o=0att  =  0 

g  g’ 

and 

Tg  =  T^,  a  =  1  at  t  *  tf 

where  tf  is  the  time  required  for  completion  of  cure.  The  degree  jf  cure 
o{t)  for  the  resin  system  can  be  determined  by  means  of  thermal  analysis.  It 
is  assumed  that  the  degree  of  cure  at  time  t  is  equal  to  the  fraction  of  heat 
released  as  measured  by  a  DSC  up  to  time  t  for  the  thermosetting  material 
under  cure.  The  total  heat  of  reaction  can  be  determined  for  the  experiments 
using  dynamic  temperature  scanning.  The  degree  of  cure  a(t)  can  then  be 
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determined  for  different  isothermal  cure  conditions  [7,8].  The  o's  thus 
determined  are  functions  of  curing  temperature  and  time. 

Equations  (1)  and  (2)  constitute  a  modified  WLF  theory  for  siit.ulating 

Chernov i scosi ty  of  thermosetting  resin.  Amotig  the  five  material  parameters,  Cj 

and  C2  can  be  determined  as  a  function  of  temperature  T  by  fitting  with 

isothermal  cure  n(a)  data,  n-r  =10^^  poises  is  given  by  Gordon  and 

9 

Simpson  [9]  for  many  polymeric  materials  at  T  =  Tg,  and  can  be  determined 
by  DSC  thermal  analysis.  In  an  isothermal  cure  case,  of  Eq.  (2)  is  equal 
to  cure  temperature.  In  the  dynamic  heating  case,  however,  is  difficult 
to  specify  exactly  and  is  the  only  adjustable  variable  in  the  present  modified 
WLF  theory. 


III.  EXPERIMENTAL  OAT/. 


For  tt.0  purpose  of  evaluating  model  predictions  of  this  modified  ULF 
theory,  experimental  data  on  Hercules  3501-6  epoxy  resin  reported  In  the 
literature  by  Lee.  Loos., and  Springer  [8]  and  Carpenter  CIO]  »ere  selected. 

These  vfill  be  sunwnarized  briefly. 

Degree  of  Cure 

Lee  et  al.  t8]  performed  thermal  analyses  on  the  epoxy  resin  system  using 
a  Differential  Scanning  Calorimeter  (OSC).  The  degree  of  cure  .  and  the  rate 
0,  degree  of  cure  do/dt  -ere  dete.-mined  from  the  results  of  the  Isothermal 
scanning  experiments.  Eight  different  temperatures  were  selected,  and  the 
results  are  reproduced  In  Figure  1.  It  is  noted  from  the  figure  that  both 
higher  values  of  the  degree  of  cure  and  rate  of  degree  ot  cure  can  be  obtained 
at  higher  cure  temperatures.  The  rate  of  degree  of  cure  Is  a  decreasing 
function  of  cure  time  t.  and  eventually  levels  off  at  a  low  rate  value  «  0.1) 
.nen  a  high  degree  of  cure  is  reached.  The  following  equations  were  found  to 
describe  these  temperature-dependent  data  accurately: 

^  =  (Ki+K2»)  o<  0.3  (3) 

da  _  If  /I  for  a  >  0.3 

where 
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\"' 


Kj  =  A,  pp  (-aEi/RT) 
K2  *  Aj  exp  (-aEj/RT) 
h  *  A3  |fexp  (-AE3/RT) 


V?lues  of  constants^ In  Eqs.  (3)  and  (4)  are  shown  In  the  Table  below. 

measured  values  of  da/dt  and  the  values  calculated  by 
Eqs.  (3)  and  (4)  are  also  shown  in  Figure  1. 


Table  1.  The  values  of  the  constants  of  Eqs.  (3)  and  (4) 
_ _ for  .‘iercules  3501-6  Resin  System  [8]  _ 

B  =  0.47 


=  2.101  X  10’  min-i 
Aj  =  -2.014  X  10’  min"* 
A3  =  1.960  X  10’  min"* 
aEi  =  8.07  X  lO**  J/mol 
aEj  =  7.78  X  10**  J/mol 
AE3  =  5.66  X  10“  J/mol 


Eqs.  (3)  and  (4)  can  be  integrated  provided  that  the  temperature  is  kept 
constant.  The  following  relationships  between  a  and  t  are  obtained: 


<2 

(1  +  ^  0.) 


b  Ln  (1-a) 


c  Ln  (1- 


(5) 


Ln  (^)  . 


(6) 
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where 


a  -  k2  (B-D/d 
2 

b  =  -  (Ki+K2)/d 
c  =  (Ki+K2B)/d 

d  =  (K1+K2)  CKiB2+(Ki-K2)  B-Kil 
and  tc  Vs  the  time  given  by  Equation  (5)  for  o*  0.3 


Chernov iscosi 


ty  n(t,T) 


CKemovtscostty  of  the  resfn  under  an  Isothennal  cure  conditton  were 
measured  by  Lee  et.  al.  tS]  on  a  Rheometrics  parallel  disc  and  plate 
mechanical  spectrometer  Model  605.  The  radius  of  the  disc  Mas  25  am 
jap  becaen  the  disc  and  the  plate  was  0.5  A  frequency  of  1.6  Hz  was 

employed.  For  a  given  cure  temperature,  the  n(t)  measured  can  be 
a(t)  obtained  from  DSC  thermal  analysis.  The  results  are  plotted  in  Figure  2 
for  three  different  Isothermal  cure  cases.  It  Is  again  noted  that  at  a  higher 
cure  temperature,  hlghc-  values  of  the  degree  of  cure  can  he  reached.  For  a 
given  «.  the  decrease  In  viscosity  Is  attrihuted  to  the  temperature  effects. 

Chemoviscosltles  of  the  3501-6  epoxy  resin  cured  under  dynamic  heating 
conditions  were  reported  by  Carpenter  [10].  Carpenter  measured  chemovisco¬ 
sltles  at  different  rate  of  Increase  In  temperatures  and  reported  the  measured 
values  as  a  function  of  cure  time.  The  experimental  data  shown  here  arc  for 
resins  cured  at  constant  heating  rates  of  l"K/m1n  (figure  5),  2.17  and  S-K/mln 
(Figure  6)  and  a  cure  cycle  (Figure  8)  similar  to  the  one  typically  found  In 

autoclave  processes  for  composite  materials. 
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The  procedure  used  for  model  evaluations  will  be  as  follows: 

1.  Model  parameters  of  Eqs.  (1)  and  (2)  will  be  determined  first  by  the 
experimental  data  of  Lee  et  al .  [8]  under  i:;othermal  cure  conditions. 
Temperature-dependent  material  constants  CjCT)  and  C2(T)  can  thus  be 
determined. 

2.  The  chemoviscosity  model  n{t,T)  established  by  isothermal  cure  data  under 
procedure  1  will  then  oe  made  to  follow  various  dynamic  heating  cure 
cycles  T(t).  The  viscosity  values  thus  calculated  will  be  compared  with 
the  experimental  data  reported  by  Carpenter  [10]. 
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IV.  HOOEL  PREDICTIONS  AND  DISCUSSION 

Isothermal  case 

Gordon  and  Simpson  [9]  pointed  out  that  glass-forming  substances 
generally  have  the  same  shear  viscosity  of  »  10^®  polses  at  their 
TgV  and  this  has  also  bee^^^^^^  polymer  systems.  By  specifying  T° 

T®  =  300®K  {27''C)  in  Eq.  (2),  the  model  (i.e.,  Eqs.  (1)  and  (2))  can  then  be 
fitted  to  the  experimental  data  n(o)  as  shown  in  Figure  2  provided  that  proper 
values  of  material  constants  C,  and  Cj  are  selected.  and  Cj  determined  for 
different  temperatures  are  tabulated  below,  and  the  model  predictions  are 
drawn  as  solid  lines  in  the  figure. 

Table  2.  Values  of  material  constants  Cj  and  C2 
in  the  modified  WLF  Theory  for  Hercules 
3501-6  resin  system  as  determined  by  the 
isothermal  cure  data. 


Cl 

CjrK) 

375 

21.1 

53.0 

400 

17.1 

33.6 

425 

18.5 

46.0 

The  model  is  capable  of  describing  the  nonlinear  relationship  between 
viscosity,  n,  and  the  degree  of  cure,  a,  as  shown  in  Figure  2,  especially  for 


the  case  of  higher  eyre  temperature.  This  is  considered  a  significant 
improvement  over  the  empirical  linear  model  used  by  Lee,  Loos,  and  Springer 
[8].  An  arrhenius  type  of  plot  is  arbitrarily  chosen  for  Cj  and  C2  as  shown 
in  Figure  3.  The  data  is  rather  scattered.  Due  to  the  (unfortunate) 
limitation  on  size  of  the  experimental  data  and  the  wide  range  in  temperatures 
covered  by  the  dynamic  heating  cure  experiments  (as  will be  discussed  later),  _ 
a  linear  least  squares  fit  was  chosen.  The  straight  lines  are  represented, 
respectively,  by: 

log  Cl  =  186.02  (1/T)  +  0.8084  (7) 

and 

log  C2  «  212.16  (1/T)  +  1.106  (8) 

It  is  noted  that  in  the  analysis  described  above,  the  glass  transition 
temperatures  Tg(t)  for  the  varying  molecular  structures  of  the  resin  cured 
at  different  tempreatures  are  chosen  as  reference  temperatures.  Ci  and  C2  are 
therefore  treated  as  "universal  constants"  in  the  present  context.  Values  of 
Cl  and  C2  reported  in  Table  2  are  close  to  the  values  of  17.44  and  51.6, 
respectively,  once  reported  by  Williams  et  al .  [3]  as  ae  universal  constants 
for  all  polymeric  materials  sytems.  Eqs.  (7)  and  (8)  will  still  be  used, 
however,  to  account  for  temperature-dependences  of  Ci  and  C2  in  chemoviscosity 
simulation  under  dynamic  heating  cure  conditions. 

The  theoretical  assumptions  as  described  above,  which  yield  close  values 
for  the  WLF  constants  Ci  and  C2,  respectively,  among  three  different 
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,sot.en.a,  caHn,  temperatures,  are  a, so  sepporteO  oi  tOe  paraUeUsm  .on, 
esper,.ota,  .ata  oOseWeO  ,n  nporo  .  .or  .  <  0.3.  Oata  „t  K,3.er  eatoes  ot 
both  0  aoO  .  are  OifficuU  to  be  measured,  and  could  be  constdered 

accurate. 


Dynamic  Heating  Case  ^ 

"  "  T„  tbe^ca«  Odating  cure  conditions,  integrations  of  E,s.  (3) 

and  (4,  are  .er,  difficult  to  perform.  If  one  assumes  that  -itbin  an,  time 
interual  at  •  t^  -  t,.  tbe  temperature  Tc  is  uniform,  and  ,s  egual  to 
,T..Ta)/3.  -Here  T.  and  are  temperatures  at  t.  and  t^  respecti.el,  as 

defined  by  the  cure  cycle,  Eqs.  (3)  and  (4)  become. 


At  =  t2  - 

=  f  (aa.TJ  - 


(9) 


where 


„4;u)-btn,..o)-cU,l.f).fof.£0.3 


and 


f  (..fa> 


for  a  >  C*- 
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the  same  as  defined  before  in  Eqs.  (5) 


The  parameters  a,  b,  c,  and  t^  are 
and  (6).  The  value  of  the  increased  degree  of  cure  02  in  Eq.  (9)  can  be 
calculated  by  the  Generalized  Newton's  Method.  Details  of  the  solution 
schemes  are  summarized  in  the  Appendix. 

Values  of  a(t)  determined  by  the  scheme  described  above  depend  heavily  on 
the  time  interval  At  chosen.  Theoretically  more  accurate  a(t)  can  be 
determined  numerically  when  a  smaller  At  is  used.  On  the  other  hand,  a 
smaller  At  could  increase  dramatically  the  numerical  round-off  errors.  A 
series  of  At  values  were  selected,  and  a(t)  was  calculated  from  Eq.  (9)  for  a 
cure  cycle  with  a  constant  heating  rate  of  5°K/min.  The  results  are  plotted 
in  Figure  4.  It  can  be  seen  that  numerical  round-off  errors  become  dominant 
when  At  <0.05  min.  It  is  also  noted  that  numerical  values  of  o(t)  calculated 
for  At  *  0.5  and  1.0  minute,  respectively,  are  almost  indistinguishable. 
Consequently,  At  =0.5  minute  was  chosen  for  all  calculations  discussed  below. 

Chernov iscosi ties  of  the  Hercules  3501-6  epoxy  resin  system  under  dynamic 
heating  cure  conditions  were  reported  by  Carpenter  [103.  Data  for  these 
constant  rates  of  heating  experiments  are  reproduced  in  Figures  5  and  6.  The 
heating  rates  are  1.0,  2.17,  and  5  ®K/min  respectively.  All  cure  cycles 
started  from  SOO^K  (27®C).  Temperature-dependent  Ci(T)  and  C2(T)  as  defined 
by  Eqs.  (7)  and  (8)  are  followed  in  the  simulation.  Values  of  the  rest  of  the 
parameters  used  in  the  model  are  tabulated  below: 


Table  3.  Values  of  model  parameters  used  in  simulating 
chemoviscosity  under  rate  of  heating  cure 
conditions 


Heating  rate 
(®K/min) 

n-j-  (poises) 

9 

T®  (“K) 

T,("K) 

1 

1013 

300 

415 

,  2.17 

10l3 

300 

430  - 

.:'-5 

1013 

300 

450 

The  model  predictions  are  shown  as  solid  curves  in  Figures  5  and  6.  It  is 
noted  that  the  theoretical  calculations  compared  favorably  with  the 
experimental  data  for  all  three  cases. 

Theoretically,  T^  is  equal  to  Tg(«)  at  T  +  «  or  a  +  1  (see  Eq.  (2)). 
Under  the  dynamic  heating  cure  conditions,  however,  T^  is  difficult  to 
specify.  Values  of  Tj;  (in  Table  3)  in  the  above  calculations  were 
selected  arbitrarily  such  that  model  predictions  are  in  reasonably  good 
agreement  with  the  experimental  data.  Nevertheless,  higher  values  of  Tc 
(Table  3)  resulting  from  the  simulation  of  higher  rates  of  heating 
experiments,  are  in  agreement  with  physical  observations  for  thermosetting 
resins. 

The  calculated  degree  of  cure  a(t)  for  all  three  cases  is  plotted  in 
Figure  7.  These  calculated  a(t)  are  based  on  the  experimental  data  for  the 
resin  under  isothermal  cure  conditions  [8],  and  are  independent  of  the 
arbitrarily  selected  values  of  discussed  above.  Values  of  da/dt  go 


through  maxima  for  all  three  cases.  The  arrow  marks  shown  in  the  figure 
denote  the  times  when  chemoviscosities  approach  infinity  during  cure.  These 
arrow  marks  are  located  approximately  at  a(t)  =  0.75  +  0.05,  and  are  near  the 
inflection  points!  of  the  curves  as  well.  Physical  significances  of  these 
observations  are  unclear  at  the  present  t^me,  and  further  research  in  the 

■I 

future  is  required.  \ 


Chernov iscosity  was  also  measured  by  Carpenter  [10]  for  the  3501-6  resin 
system  under  a  typical  cure  cycle  used  in  autoclave  processing  for  composite 
materials.  The  experimental  data  is  reproduced  in  Figure  8.  Values  of  model 
parameters  used  are  the  same  as  those  included  in  Table  3  and  Eqs.  (7)  and 
(8),  except  that  Tc  =  375“K.  The  theoretical  predictions,  shown  as  solid 
curves,  are  seen  in  reasonably  good  agreement  with  the  experimental  data.  The 
model,  however,  fails  to  describe  the  chemoviscosity  "shoot-up"  near  the  end 
of  the  cure  cycle.  The  reasons  are  unclear  at  the  present  time.  It  could  be 
due  to  the  inadequate  definitions  of  Ci(T)  and  C2(T)  by  Eqs.  (7)  and  (8),  or 
due  to  different  reaction  mechanisms  appeared  in  the  high  curing  temperature 
regime.  Calculated  values  of  «(t)  are  plotted  in  Figure  9.  The  arrow  mark 
shown  denotes  the  time  when  '.hemov iscosity  "shoot-up"  occurs.  Features  of  the 
o(t)  curve  discussed  above  for  isothermal  cure  cases  are  again  noticeable  in 

this  figure. 


Temperature  ranges  covered  by  Carpenter  [10]  in  dynamic  heating  cure 
experiments,  as  shown  by  dotted  lines  in  Figure  3,  far  exceed  the  ranges 
covered  by  Lee  et  al .  [8]  in  isothermal  cure  experiments  (shown  by  solid  lines 
in  the  figure).  Accurate  definitions  of  the  crucial  temperature-dependences 


#.Ci{T)  and  C2(T)  needed  for  the  model  thus  suffer  from  serious  data 
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extrapolattoa  errors.  The  results  of  the  preseot  Investtgatloos,  although 
eocouraglhg,  cah  ooly  de»hstr.te  the  feaslhilit,  of  the  approach  of  nodifyiog 
the  «LF  theory  io  checoviscosit,  modellog  for  theruiosettlog  reslo.  Detailed 
studies  such  as  the  precise  definitions  of  Tg(a)  (Tg*  (21)»  dhd  the 
functions  of  C,(T)  and  Cj(T)  over  the  full  temperature  range  by  experimental 
measurements  for  the  resin  system  under  Investigation  have  to  be  pursued  in 

the  future. 
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V.  CONCLUSIONS 

A  new  analytical  model,  based  upon  modifications  of  the  WLF  theory,  has 
been  formulated.  The  model  assumes  that  the  glass  transition  temperature 
Tg(t)  of  the  resin  system  under  cure  is  linearly  related  to  the  degree  of 
cure  a(t)  determined  by  DSC  thermal  analysis.  Theoretical  predictions  or^ 
Chernov iscosity  compare  favorably  with  the  experimental  data  reported  in  the 
literature  for  both  isothermal  and  dynamic  heating  cure  conditions. 
Considering  that  data  are  available  only  in  the  relatively  narrow  temperature 
range  used  in  isothermal  curing  experiments  by  Lee,  Loos,  and  Springer  [8], 
the  agreement  between  theoretical  predictions  and  experimental  data  reported 
by  Carpenter  [10]  for  dynamic  heating  cure  cycles  are  particularly 
encouraging.  Due  to  the  limitations  in  size  of  the  data  available,  the 
successful  analysis  reported  here  has  established  the  feasibility  of  the 
proposed  approach  which  simulates  chemoviscosity  for  thermosetting  resin  by 
modifying  the  widely  used  WLF  theory  in  polymer  rheology. 
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Appendix  -  Generalized  Newton's  Method  in  Solving  02  from  Eq.  (9) 

Eq.  (9)  can  be  rearranged  to  become: 

9(02)  =  f(a2»T^)  -  f  {ai,T^)  -  At^  =  0  (a) 


and.  ^ 

g'Caa)  =  f’(a2.Tc^  (b) 

where  g'(o2)  =  dg/da2.  and  f (a2.Tc)  =  df/da2»  Gentralized  Newton's 
Method  sets  up  the  following  recu.sion  formula  to  solve  for  02  in  Eq.  (a): 

~  ®2  ~  w.g(o2)/g ' (012) »  ~  0»  if  2...  (c) 

provded  that  g'(a2)/*0;  where  w  is  called  the  over-relaxation  factor.  A  good 
choice  for  w  can  be  determined,  in  general,  cnly  by  experimentation.  A  choice 
of  w  different  from  unity  often  can  increase  the  convergence  rate  of  the 
routine  appreciably.  In  this  study,  w  =  0.1  was  used  for  all  calculations. 
Iterations  are  continued  until  -  a"|<  lu®. 
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Figure  1.  Rate  of  degree  of  cure  vs.  degree  of  cure 
reported  by  Lee,  Loos,  and  Springer 
Hercules  3Sni-6  resin  system  under 
isothermal  cure  cycles.  Solid  lines  were 
calculated  using  Eqs.  (3)  and  (4). 


a,  degree  of  cure 


Figure  2.  Chemoviscosity  as  a  function  of 

degree  of  cure  for  Hercules  3501-6 
resin  system  under  isothermal  cure 
cycles.  Symbols  represent 
experimental  data  reported  by  Lee, 
Loos,  and  Springer  [8].  Solid 
lines  were  calculated  by  the  model , 
Eqs.  (1)  and  (2). 


3  Temperature-dependences  of  the  WLF 

constants  Cl  and  Cj  determined  by  the 
isothermal  cure  data  of  Lee  et  al .  I.o.l 
for  the  Hercules  3501-6  resin  system. 


Figure  4.  Simulated  degree  of  cure  data  for  the  Hercules 
3501-6  resin  system  cured  under  constant  rate 
of  heating  condition,  5“K/min.  Various  time 
intervals  At  were  used  in  solving  Eq.  (9). 
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Figure  6.  Chernov iscosity  of  the  Hercules  3501-6  resin 
system  cured  under  constant  rates  of  heating 
conditions  of  2.17  and  B^K/min,  respectively. 
Circles  denote  the  experimental  data  by 
Carpenter  [10].  Solid  curves  denote  the  model 
predictions  of  Eqs.  (1)  and  (2). 


t,  min 


Figure  8.  Chemoviscosity  of  the  Hercules 

3501-6  resin  system  cured  under  the 
cure  cycle  shown.  Circles  denote 
the  experimental  data  by  Carpenter 
[10].  Solid  curve  was  calculated  by 
the  model,  Eqs.  (1)  and  (2). 
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